Levy-flight intermixing: anomalous nanoscale diffusion in Pt/Ti 
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Probing the anomalous nanoscale intermixing using molecular dynamics (MD) simulations in 
Pt/Ti bilayer we reveal the superdiffusive nature ol interfacial atomic transport. It is shown that 
the Pt atoms undergo anomalous atomic transport across the anisotropic interface of Pt/Ti with 
suprisingly high rates which can be characterized as Levy flights. Levy flight is not a unique 
phenomenon in nature, however, no such events have been reported yet for bulk interdiffusion. In 
particular, the low-energy (0.5 keV) ion-sputtering induced transient enhanced intermixing has been 
studied by MD simulations. Ab initio density functional calculations have been used to check and 
reparametrize the employed heteronuclear interatomic potential. The Levy-intermixing behavior 
explains the high diffusity tail in the concentration profile obtained by Auger electron spectroscopy 
depth profiling (AES-DP) analysis in Pt/Ti bilayer (reported in ref.: P. Siile, et al., J. Appl. Phys., 
101, 043502 (2007)). 

PACS numbers: 05.40.Fb, 68.35. Fx, 02.50.Ey, 66.30.Jt, 61.80.Jh, 66.30.-h, 



Anomalous diffusion has long been known since the 
classical Richardson's study on turbulent diffusion El- 
Since than anomalously large diffusion rates have been 
observed in various systems under different conditions 
@> B B IE 0> HI • The common feature of the observed 
enhanced diffusion rates are the nonlinear growth of the 
mean squared displacement in the course of time leading 
to a power-law time scaling pattern Hence we are 
faced with non-Brownian motion of particles with non- 
Gaussian diffusion rates giving rise to anomalous behav- 
iors which are termed as strange kinetics [9(. Superdif- 
fusion takes place with Levy flight characteristics (LF, 
Levy jump diffusion with broad spatial jump distribu- 
tions) above the Brownian motion regime [2]. In patic- 
ular, LF behaviour has been found in turbulent mixing 
in fluids (review in refs. @, B 0) or m elongated mi- 
celles [6|, for photon trajectories in incoherent atomic 
radiation trapping [Ic] ] and it became clear that Levy 
processes can be important in the analysis of chaotic dy- 
namics 0]. The research of LF behaviors covers many 
disciplines of natural and social sciences hence this topic 
is of high interest @, IE 0, IE 0, El • Our primary purpose 
is to present experimental and computational evidences 
that the nanoscale transient enhanced intermixing and 
diffusion, a common phenomenon in many nanostruc- 
tured materials, could also be fitted in the category of 
Levy flight superdiffusion. 

There are a growing number of evidences emerged in 
the last decades that anomalous nanoscale broadening 
of interfaces or high diffusity tail in the impurity con- 
centration profile occur during ion-irradiation or sput- 
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14 . [lil Hit ], sputter deposition and thin 



tcring 

film growth |17l . [18( . The anomalous nanoscale bulk dif- 
fusional effects could be due to still not clearly estab- 
lished atomistic accelerative effects leading to anoma- 



lously fast and possibly athermal atomic transport (long 
ran ge d iffusion) at interfaces or during impurity diffusion 

mM 0, EE EE El HI . 



In the present work, using computer atomistic simu- 
lations, the recently observed enhanced intermixing in 
a nanoscale bilayer film/substrate system (Pt/Ti) re- 
ported in ref. |l6[ is attempted to interpret as a Levy- 
superdiffusive atomic transport process in the bulk. Su- 
perdiffusion in the nanoscale has only been reported until 
now on solid surfaces 2l|, l22| and no reports have been 
found for bulk superdiffusion in the solid state of matter. 

Classical molecular dynamics simulations have been 
used to simulate the ion-sputtering (repeated low-energy 
ion bombardments of the film) induced atomic intermix- 
ing at the Pt/Ti interface using the PARC AS code 
Here we only shortly summarize the most important as- 
pects. A variable timestep and the Berendsen tempera- 
ture control is used to maintain the thermal equilibrium 
of the entire system at 300 K (see refs. at 16]). The 
detailed description of other technical aspects of the MD 
simulations are given in refs. [23| and details specific to 
the current system in recent communications 

EE Hi Hi- 

Our primary purpose is to simulate the conditions oc- 
cur during ion-sputtering [1 61 ] and Auger electron spec- 
troscopy depth profiling analysis (AES-DP) EE using 
molecular dynamics simulations. Following our previous 
work EH we ion bombard the film of the bilayers Pt/Ti 
and Ti/Pt (film/substrate systems) with 0.5 keV Ar + 
ions repeatedly (consequtively) with a time interval of 
10-20 ps between each of the ion-impacts at 300 K which 
we find sufficiently long time for most of the structural 
relaxations and the termination of atomic mixing, such 
as sputtering induced intermixing (IM) [251 ]. Since we 
focus on the occurrence of transient intermixing atomic 
transport processes, the relaxation time of 10 — 20 ps 
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TABLE I: The semiempirical parameters used in the tight 
binding interatomic potential for Pt/Ti [ltj |26| 
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Ti 


1.416 


1.643 


0.074 


11.418 


2.95 


Pt 


2.695 


4.004 


0.298 


10.612 


2.78 


Ti-Pt 


4.2 


2.822 


0.149 


11.015 


2.87 





"The parameters for Ti and Pt have been given by Cleri and 
Rosato |2tSll . The parameters of the crosspotential have been ob- 
tained by fitting the interpolated crosspotential given in ref. [l6l | 
to ab initio diatomic calculations. The preexponential parameter 
£ has been fitted to the ab initio curve (see also in ref. jl9j|'l. The 
notations used for the parameters are the same as given in refs. 

ee mm. 

should be appropriate for getting adequate information 
on transient enhanced intermixing. Pair potentials have 
been used for the interaction of the Ar + ions with the 
metal atoms derived using ab initio density functional 
calculations. The initial velocity direction of the impact- 
ing ion was 10 degrees with respect to the surface of the 
film (grazing angle of incidence) to avoid channeling di- 
rections and to simulate the conditions applied during 
ion-sputtering (lfj| . The impact positions have been ran- 
domly varied on the surface of the film/substrate system 
and the azimuth angle (f> (the direction of the ion-beam) . 
In order to simulate ion-sputtering a large number of ion 
irradiation are applied using script governed simulations 
conducted subsequently together with analyzing the his- 
tory files (movie files) in each irradiation steps. The im- 
pact positions of the 400 ions are randomly distributed 
over a 20 x 20 A 2 area on the surface. 

We used a tight-binding many body potential on the 
basis of the second moment approximation (TB-SMA) to 
the density of states [26| , to describe interatomic inter- 
actions. In ref. [lfl it has been shown that the TB- 
SMA potential gives the reasonable description of IM 
in Pt/Ti and gives interfacial broadening comparable 
with AES-DP measurements. The crosspotential energy 
has been calculated for the Ti-Pt dimer u sing ab initio 
local spin density functional calculations [27j together 
with quadratic convergence self-consistent field (SCF) 
method. The G03 code is well suited for molecular calcu- 
lations, hence it can be used for checking pair-potentials. 
We used the Perdew-Burke-Ernzerhof (PBE) gradient 
corrected exchange-correlation potential J28| . We find 
that the interpolated TB-SMA potential [16| nearly per- 
fectly matches the ab initio one hence we are convinced 
that the TB-SMA model accurately describes the het- 
cronuclear interaction in the Ti-Pt dimer. The employed 
parameter set is given in Table 1. 

The crossectional computer animations of simulated 
ion-sputtering can be seen in our web page [29( . The car- 
toons of the simulation cells (crossectional slabs as a 3D 
view) can be seen in ref. [l6| which show the strong inter- 
mixing at the interface in Pt /Ti and a much weaker mix- 



ing in Ti/Pt. In Fig [I] the evolution of the mean-square 
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FIG. 1: The simulated mean square of intermixing atomic 
displacements (R 2 ) z (MSD, variance in A /particle) in Pt/Ti 
and Ti/Pt as a function of the ion-fiuence (the number of ion 
bombardments, Ni on ) obtained during the ion-sputtering of 
these bilayers at 500 eV ion energy (results are shown up to 
~ 400 ions). Results from two independent simulations are 
shown for Pt/Ti. Fits to the variance of various curves are 
also given with dotted lines which is obtained as {R 2 )z ~ N a . 
The obtained values of a are also shown. 

of atomic displacements (MSD) (the variance of displace- 

ments) (R 2 (t)) z = £f •» ^— £f [rj (t) - r ; (t = 

0)] 2 , of all intermixing atoms obtained by molecular dy- 
namics simulations, where (iV (£) is the position vector of 
atom 'i' at time t, N 3 atom {t) is the total number of inter- 
mixing atoms in the jth irradiation step included in the 
sum), can be followed as a function of the ion fluence (the 
number of ions Ni on ). Lateral components (x,y) are ex- 
cluded from (R 2 ) z and only contributions from IM atomic 
displacements perpendicular to the layers are included (z 
or depth components). We follow during simulations the 
time evolution of (R 2 ) z which reflects the atomic migra- 
tions through the interface (no other atomic transport 
processes are included). 

In Fig. Q]we present (R 2 ) z as a function of the number 
of ion impacts N ion (ion-number fluence). (R 2 ) z {N ion ) 
corresponds to the final value of (R 2 ) z obtained during 
the N ion t\i simulation. The final relaxed structure of 
the simulation of the (iVj on — l)th ion-bombardment is 
used as the input structure for the A^ ort th ion-irradiation. 
The asymmetry of mixing can clearly be seen when 
{R 2 ) z (Ni on ) and the depth profiles given in ref. [HI are 
compared in Ti/Pt and in Pt/Ti. The computer anima- 
tions of the simulations [2^] together with the plotted 
broadening values at the interface in ref. [li} also re- 
veal the stronger IM in Pt/Ti. Moreover we find the 
strong divergence of (R 2 ) z from linear scaling for Pt/Ti 
and a much weaker nonlinear scaling has been found for 
Ti/Pt. As it has already been shown in ref. [l6| AES- 
DP found a relatively weak IM in Ti/Pt (the interface 



FIG. 2: The normal to the surface crossectional view of 
collisional displacement cascades with atomic trajectories 
(crossectional slab cut in the middle of the simulation cell) 
in Pt/Ti and in Ti/Pt. The right and left panels correspond 
to Pt/Ti and to Ti/Pt, respectively. The positions of the en- 
ergetic particles of which the kinetic energy is larger than ~ 1 
eV are collected during typical single ion impact events at 
500 eV ion energy. The ions have been initialized 10 A above 
the interface. The vertical axis corresponds to the depth po- 
sition given in A. The position at z — is the depth position 
of the interface with a horizontal dashed line and the initial 
positions of the atoms (before irradiation) are also shown. 



FIG. 3: The normal to the surface crossectional view of tra- 
jectories of Levy flights obtained during 50 consequtive ion- 
impact events in Pt/Ti using 0.5 keV ion energy for Pt (left 
panel). The atomic trajectories of hyperthermal Ti atoms 
are also shown collected during the same irradiation events 
(panel on the right). The atomic layers above the interface 
with depth position z > are Pt layers. The positions of en- 
ergetic atoms have been plotted which have larger than ~ 2 
eV kinetic energy. The ions have been initialized 10 A above 
the interface. The interface is at z — depth position (the in- 
crements are in A). The initial positions of the atoms (before 
irradiation) are also shown. 



broadening a ss 20 A) while an unusually high IM occurs 
in the Pt/Ti bilayer (a ~ 70 A). In order to clarify the 
mechanism of intermixing and to understand how much 
the nanoscale interfacial mass-anisotropy influences IM, 
simulations have been carried out with atomic mass ra- 
tio S — mpt/rriTi (where mpt and m-r, are the atomic 
masses) is artificially set to S ~ 1 (mass-isotropy). We 
find that the magnitude of IM is strongly sensitive to 6 at 
mass-anisotropic, called (^-interfaces. The corresponding 
animation can be also be seen at [29[ . We reach the con- 
clusion that the mass-effect is robust and the magnitude 
of IM is weakened significantly in artificially isotropic 
Pt /Ti in agreement with our earlier finding (2Jj [25] . The 
huge difference in IM between Pt/Ti and Ti/Pt can be 
understood as the effect of ^-inversion on IM. Actually 
the system undergoes the transition in the asymptotics 
of (R 2 ) oc t 135 — > (R 2 ) oc t 1 - (from nonlinear dynamics 
to linear) when the mass anisotropy is inverted (the film 
and the substrate is interchanged). 

To further test mass-effect on IM, we carried out simu- 
lations for the Pt/Ti system in which the atomic masses 
have been interchanged (Ti possesses the atomic mass of 
Pt and vice versa) setting in an artificial mass ratio (the 
inverse of the normal one) while keeping all the other pa- 
rameters unchanged. Note, that we keep the interatomic 
potentials, only the atomic mass ratio is inverted. We 
find that this artificial setup of atomic masses results in 
the suppression of IM in Pt/Ti with inverted 6. These 
findings together with our AES measurements (with the 
long-range tail shown in ref. [IB]) confirms our recent 
results reported for various bilayers in which a strong 
correlation has been obtained between the experimental 



and simulated mixing efficiencies and mass anisotropy in 
various metallic bilayers [2~i ]. 

In Fig. [T]for Pt/Ti it can clearly be seen that (R 2 ) z 
(MSD) grows nonlinearily with Ni on . The horizontal 
axis is proportional to the time of ion-sputtering, hence 
(R 2 ) z °c t a with a > 1. [2j. During thermally activated 
interdiffusion {R 2 ) scales linearily with time hence our 
system violates this temporal behavior. The (R 2 ) z oc t a , 
scaling, where a > 1, used to be considered as the signa- 
ture of anomalous diffusion (superdiffusion) in the liter- 
ature @, 0, HI- In our case we find that a ~ 1.35 ± 0.1 
exponent fits the simulated curve for Pt/Ti. We find 
for Ti/Pt a linear behavior (a m 1.0 ± 0.05) hence no 
superdiffusive features can be identified in this system. 

We would like to show that it might also be the case 
that transient IM takes place in Pt/Ti which resembles 
in many respect the Levy superdiffusive atomic trans- 
port processes known on solid surfaces Q . Superdiffu- 
sion occur when the particle trajectories exhibit long dis- 
placements (flights) termed Levy flights Q , characterized 
by power-like asymptotical distribution of the displace- 
ments of trajectories (and with heavy-tailed distribution) 
0. For these ballistic processes {R 2 )(t) is divergent with 
time. Normal diffusion is characterized by the linear time 
evolution of (R 2 ) in the long-time limit and with Gaus- 
sian distribution [3j. 

The obtained results are in agreement with the findings 
presented in ref. (l6| . However, in that paper it has not 
been realized that the fingerprint of superdiffusive fea- 
ture of IM is detected by AES as a long-range diffusity 
tail in the concentration profile at the Pt/Ti interface 
in the Pt/Ti bilayer. No such tail occurs in the concen- 
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tration profile of Ti/Pt as shown in ref. Hence we 

conclude that the experimental fingerprint of Levy jumps 
could also be detected as a high diffusity tail in the AES 
concentration profile of many other anisotropic materials. 

The trajectories corresponding to Levy flights can be 
seen in the left panel in Fig. O In the panel of Pt/Ti in 
Fig [2] we can see the ballistic trajectories of intermixing 
hyperthermal Pt atoms which exhibit Levy flights with 
long trajectories through the interface. We find 8 events 
out of 10 which exhibit Levy flights. No such trajectories 
can be seen in the panel of Ti/Pt in Fig. [5] and in 
the other dozens of events generated (not shown). The 
ions have been initialized during these single ion impact 
events 10 A above the interface in order to place the range 
of the ions directly in the depth position of the interface. 
This way of direct deposition of the ion energy at the 
interface does not influence the main physics what we find 
with simulations when ions have been initialized at the 
surface. This is because the first few tens of ions simply 
sputter remove the top layers of the film and intermixing 
develops when the range of ions approaches the interface. 
Fig. [3] depict us what we see in Fig. [TJ that Levy flights 
cause the nonlinear time scaling of (R 2 ) z - The initial 
kinetic energy of few of these particles can reach the few 
tens of eV and which hyperthermal Pt atoms exhibit long 
trajectories. 

In the left panel of Fig. [3] we show all the flight tra- 
jectories of Pt atoms obtained during a simulation with 
50 repeated ion impacts with 10 ps time delay between 
each of the events in Pt/Ti. This figure clearly depicts us 
that Levy jumps of energetic Pt atoms take place which 
boost intermixing. Such kind of an interfacial broaden- 
ing can also be called Levy interdiffusion. In the right 
panel of Fig. [3J the trajectories of energetic Ti atoms 
are also shown obtained during the same simulation. No 
Levy jumps can be seen for Ti atoms. These atoms not 
even move across the interface except during few events. 
In Ti/Pt no LF behavior can be found both for Pt and 
Ti atoms. Concerning the atomic mobility of the hyper- 
thermal Pt particles, some of these particles can have a 
huge initial kinetic energy exceeding 100 eV in few cases. 
This could be due to the accelerative effect of head on 
collision of the ion or recoils with few of the Pt atoms 
at the 5-interface. The Pt atoms intermix preferentially 
both in Ti/Pt and in Pt/Ti [25j. These rapidly migrating 
particles slow down to few eV within 0.5 ps and which 
energy regime persists up to few ps. Hence a long range 
tail exists for the time distribution of the atomic velocity 
during Levy jumps in Pt/Ti which is also a characteristic 
feature of non-Brownian dynamics 0, Q • 

In conclusion, we reveal that atomic intermixing might 
take place via Levy flights in the Pt/Ti film/substrate 
bilayer upon low-energy ion bombardments. This mech- 
anism could be valid for many other interdiffusion pro- 
cesses in which mass anisotropic interface is present. 

This work is supported by the OTKA grants F037710 



and K-68312 from the Hungarian Academy of Sciences. 
The work has been performed partly under the project 
HPC-EUROPA (RII3-CT-2003-506079) using the super- 
computing facility at CINECA in Bologna and at NIIF 
(Budapest). 



[1] L. F. Richardson, Proc. Roy. Soc. (London) Ser. A110, 
709 (1926). 

[2] R. Metzler, J. Klafter, Phys. Rep., 339, 1. (2000), Levy 
Flights and Related Topics in Physics, ed. by M. F. 
Shlesinger, G. M. Zaslavsky, and U. Frisch (Springer, 
Berlin, 1995). 

[3] G. M. Zaslavsky, in Hamiltonian Chaos & Fractional Dy- 
namics, p. 229 (Oxford, 2005). 

[4] R. van Gastel, et al, Nature, 408, 665 (2000), Phys. Rev. 
Lett. 86, 1562 (2001). 

[5] A. M. Edwards, et al, Nature, 449, 1044 (2007). 

[6] J. P. Bouchaud, et al, J. Phys. II. (France) 1, J. P. 
Bouchaud, A. Georges, Phys. Rep. 95, 127 (1990). 

[7] B. I. Shraiman, and E. D. Siggia, Nature (London) 405, 
639 (2000). 

[8] J. Tsai, et al, Nano Lett. 8, 425. (2008). 
[9] M. F. Shlesinger, G. M. Zaslavsky, and J. Klafter, Nature 
363, 31 (1993). 

[10] E. Pereira, J. M. G. Martinho, M. N. Berberan-Santos, 
Phys. Rev. Lett. 93, 120201-1 (2004). 

[11] J. Klafter, M. F. Shlesinger, and G. Zumofen, Phys. To- 
day 49, 33 (1996). 

[12] G. Abrasonis, W. Moller, X. X. Ma, Phys. Rev. Lett. 96, 
065901-1, (2006). 

[13] N. Bibic, et al, Appl. Phys. Lett. 90, 051901 (2007). 

[14] J. Cardenas, et al, Appl. Phys. Lett., 73, 3088 (1998). 

[15] W. Nieveen, et al, Appl. Surf. Sci., 231-232, 13067 
(1989). 

[16] P. Side, et al, J. Appl. Phys., 101, 043502 (2007). 

[17] J. D. R. Buchanan, et al, Phys. Rev. B66, 104427 

(2002) . 

[18] P. Sule, J. Chem. Phys., (2008), 

arXiv:condmat/0802.0347. 
[19] P. Sule, submitted for publication, 

arXiv:condmat/0802.0427. 
[20] D. K. Avasthi, et al, Phys. Rev. B68, 153106 (2003). 
[21] W.D. Luedtke and U. Landman, Phys. Rev. Lett. 82, 

3835 (1999). 

[22] J. G. Amar, F. Family, D. C. Hughes, Phys. Rev. B58, 
7130 (1998), O. V. Bychuk, B. O'Shanghnessy, Phys. 
Rev. Lett. 74, 1795 (1995). 

[23] K. Nordlund, et al, Phys. Rev., B57 7556. (1998). 

[24] P. Sule, M. Menyhard, Phys. Rev., B71, 113413 (2005). 

[25] P. Stile, M. Menyhard, K. Nordlund, Nucl Instr. and 
Meth. in Phys. Res., B226, 517 (2004), B211, 524 

(2003) . 

[26] V. Rosato, M. Guillope, B. Legrand, Phil. Mag. A59, 321 
(1989), F. Cleri and V. Rosato, Phys. Rev. B48, 22-33 
(1993). 

[27] M. J. Frisch, G. W. Trucks, H. B. Schlegel, et al, 
Gaussian, Inc., Pittsburgh PA, (2003), see also at: 
http:/ /www. gaussiaii.com 

[28] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[29] http://www.mfa.kfki.hu/~sule/animations/ptti.html 



